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Abstract 



We use data collected with the CLEO II detector to perform a high-statistics 
measurement of the resonant substructure in K^tt^tt^ decays. We 

find the Dalitz Plot is well represented by a combination of seven quasi- 
two-body decay channels (K*°tt°, K' p, K*-tt+, iro(1430)-7r+, Ko(1430)°7r°, 
K~ p~^(1700), and i^*(1680)~7r~''), plus a small non-resonant component. We 
see no evidence of a scalar k — > K~7r~^ resonance in the mass range recently 
reported by other groups. Using the amplitudes and phases from this analysis, 
we calculate an integrated CP asymmetry of —0.031 it 0.086. 
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I. INTRODUCTION 



A clearer understanding of final state interactions in exclusive weak decays is an important 
ingredient for our ability to model decay rates as well as for our understanding of interesting phe- 
nomena such as mixing Q. There are several theoretical methods used to understand the 
dynamics of two body charmed meson decays with experimental measurements as input. Unfortu- 
nately, final-state interactions are often not well understood, and may not be included properly in 
many models. These long-distance strong interaction effects can cause significant changes in decay 
rates for specific final states, and can cause shifts in the phases of the decay amplitudes. 

Three-body decays provide a rich laboratory in which to study the interference between in- 
termediate state resonances, and provide a direct probe of the final state interactions in certain 
decays. When a particle decays into three or more daughters, intermediate resonances dominate 
the decay rate. These resonances will cause a non-uniform distribution of events in phase space 
when analyzed using a "Dalitz Plot" technique Q. Since all events of a particular decay mode 
have the same final state, multiple resonances at the same location in phase space will interfere. 
This provides the opportunity to experimentally measure both the amplitudes and phases of the 
intermediate decay channels, which in turn allows us to deduce their relative branching fractions. 
These phase differences can even allow details about very broad resonances to be extracted by 
observing their interference with other intermediate states. 

This paper describes a study of the underlying structure in K^tt^tt^ decays Four 

previous groups have made similar measurements, each with less than 1000 signal events ||8-11|. In 
addition to significantly increasing our statistical power, the large CLEG II dataset on which this 
analysis is based also permits us to tighten analysis requirements to drastically reduce the effect of 
backgrounds. Coupled with the superb resolution of the CLEO-II detector, this has allowed us to 
extract significantly more information about this decay than was possible in the past. 



II. THEORETICAL MODELS 

Since we are studying the decay of a spin-zero particle to three daughters, only two degrees- 
of-freedom are required to completely describe the kinematics. To see this, consider the decay 
in the rest frame. The four-momenta of the three final state particles correspond to twelve 
unknowns. We have one constraint for each known mass and four additional constraints from the 
conservation of momentum and energy in the decay. Finally, since the three degrees of freedom 
describing the spatial orientation of the decay are irrelevant (the having spin zero) only two 
independent degrees of freedom remain. 

There are three invariant masses that can be formed by considering all possible pairs of final 
state particles: M|,_^+, M'^_^o and M^+^o- Only two of these are independent, however, since 
energy and momentum conservation results in the additional constraint 

Mlo + M|_ + + M^o = M|_^+ + M|_^o + M^o- (1) 

Choosing two of the above three invariant masses as dynamic variables has two compelling 
advantages: i) Their relativistic invariance means the Lorentz frame in which they are evaluated 
is irrelevant; ii) As the expression for the partial width in Equation |2| indicates, we expect that a 
scatter plot of events in the vs -/Vf23 plane (known as a Dalitz Plot) will be uniformly distributed 
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FIG. 1. A Representation of the three body decay of — > ABC through an AB resonance. 
The spin sum is performed to obtain the angular dependence of the decay. 

if phase-space alone determines the decay dynamics. This allows the decay fraction at each point 
to be readily correlated with the decay matrix element M without additional corrections: 

In this analysis we choose M^_^+ and M^4_^o as our two Dalitz Plot variables, and must next 
construct the relevant decay amplitudes in terms of these. The dynamics can be understood with 
the aid of Figure |l[ We first consider the decay of the meson into particle C plus an AB 
resonance, followed by the decay of the AB resonance into particles A and B. To properly describe 
the structure of this decay using our Dalitz Plot variables, we need to obtain the angular dependence 
of the decay products. Each vertex in Figure || contains a spin factor £\ which depends on the type 
of the decay: scalar, vector, tensor, etc. The matrix element for a vector decay is 



' - Mis - iMrV ab' 



M = FoiPDO + ^c)m ..9 ..2 ..T. (Pa - PB)uFr (3) 



where P denotes 4-momentum, and is the mass of the resonance. In general the form factors 
at each vertex, F£) and F^, are unknown functions, however in practice they are either set to a 
constant value or to the Blatt-Weisskopf penetration factors [12|. 



The spin-sum in the numerator of Equation |3| is evaluated to give 

T.^Te\ = -r + ^^^ (4) 

A ^^^AB 

and the "mass dependent width" Tab is a function of the AB invariant mass M^b-, the momentum 
of either daughter in the AB rest frame pab, the momentum of either daughter in the resonance 



rest frame pr, the spin of the resonance J, the width of the resonance F.^, and is expressed as [13|: 

2J+1 
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We relax the transversality requirement on the vector resonance in Equation |^ and divide by 
instead of M^^. This substitution gives rise to a small spin zero component when the vector 

resonance is off mass-shell, a behavior which is observed to occur with the W boson and which 

should also be expected in the resonance behavior we are studying here. 

Inserting this expression for the spin-sum into Equation ^ and summing over the repeated indices 

gives the Lorentz invariant expression for the matrix element of a vector particle as a function of 

position in the Dalitz Plot: 



Ai{ABC\r) = FoFr 



2 {Ml-M}^(Ml-Ml) 



M2 - M\j^ - iMrVAB 



(6) 



The procedure for calculating the vector matrix element is generalizable to intermediate particles 
having other spin. For example, we can easily find the amplitude for a spin zero resonance to be 



Ao{ABC\r) = FoFr 



1 



M^-Mls-iMrTAB' 



(7) 



The procedure for higher spin resonances involves a bit more algebra. For example, the spin 
two case starts with 



A2iABC\r) = Fd{Pd + Pc)i,{Pd + Pc\ 



M2 - Mis - iMrVAB 



{PA-PB)a{PA-PB)pFr. (8) 



In this case the spin sum has been previously calculated by Pilkuhn |13| to be 



E 



2 V 3 



(9) 



where 



rp/lU 



(10) 



When this expression is inserted into Equation ^ and the implied sums performed we find the 
final form of the tensor matrix element: 



A2{ABC\r) 



FnFr 



M2 - - iFABMr 



M'ab 



2Ml 



2M1; + 



(M 



M2^2 



M2 




BC 



M 



AC 



+ 



{Ml-Ml)(Ml-Ml) 



M2 



MIb 



2M\ 



2M| + 



Ml?' 



M2 



(11) 



Next we return to the form factors and -Fj., which attempt to model the underlying quark 
structure of the meson and the intermediate resonances. We use the Blatt-Weisskopf penetration 
factors shown in Table |. These have one free parameter, R, which is the "radius" of the meson, 
and are dependent on the momentum P of the decay particles in the parent rest frame. In all cases, 
we normalize the form factor to have unit value at the nominal meson mass. The fits display very 
little sensitivity to the meson radii; good fits are obtained when these values vary between GeV~^ 



and 10 GeV ^ for the and between GeV ^ and 3 GeV ^ for the intermediate resonances. To 
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TABLE I. Blatt-Weisskopf Penetration Form Factors, pr is the momentum of either daughter 
in the meson rest frame. pAB is the momentum of either daughter in the candidate rest frame 
(same as pr except the parent mass used is the two-track invariant mass of the candidate rather 
than the mass of the meson). R is the meson radial parameter. 



Spin 


Form Factor 





1 


1 




2 


-y/9+3K2p2_,_^4p4 







be consistent with other experiments Q we have chosen the to have R = 5 GeV~^ and the 
intermediate resonances all to have R = 1.5 GeV~^ 

Before continuing, we must specify our phase conventions for the intermediate resonances. We 
can explicitly see the importance of specifying the ordering of particles in the decay by examining 
Equation |6[ If we were to switch the labels A and B we would generate an overall minus sign 
causing the phase to change by 180°. In an attempt to be consistent with previous results we have 
chosen the phases in the same way as the E687 collaboration |^ since they are the only group to 
have explicitly published their choice of phases and matrix elements. 

Now that we know the form of the intermediate resonance amplitudes, and have chosen a phase 
convention that will allow us to compare our results with previous measurements, we can write 
down an expression for the overall matrix element of the decay. Guided by the results of previous 



measurements [|-1C], we begin with only three vector resonances p{770)^ , K* and K*^ [15| as 
well as a flat non-resonant (nr) component: 

+ape^'^Mi(7r+7r°i^-|p+) 

+aK*-e''^'<'- Ai {K-Tr\+\K*- ) , (12) 

where the ai and are the amplitude and relative phase of the i'th component respectively. The 
overall normalization is arbitrary, and is chosen to be 

\M\'^dVV = 1 (13) 

where dT>V indicates that the integral is performed over the Dalitz Plot. This is equivalent to 
saying that we are sensitive only to relative phases and amplitudes, which in turn means that we 
are free to fix one phase and one amplitude in Equation [l^. To minimize correlated errors on 
the phases and amplitudes we choose the largest mode, K~ p, to have a fixed zero phase and an 
amplitude of one. 
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Since the choice of normahzation, phase convention, and amphtude formahsm may not always 
be identical for different experiments, fit fractions are reported instead of amplitudes to allow for 
more meaningful comparisons between results. The fit fraction is defined as the integral of a single 
component divided by the coherent sum of all components: 



/ are''f'-A{ABC\r] 

Fit Fraction = ■ 



/ j:,aje^'l'^A{ABC\j: 



dVV 

-2 • (14) 

dVV 



The sum of the fit fractions for all components of a fit will in general not be one because of 
interference. 

One must also consider the statistical errors on the fit fractions. We have chosen to use the 
full covariance matrix from the fits to determine the errors on fit fractions so that the assigned 
errors will properly include the correlated components of the errors on the amplitudes and phases. 
After each fit, the covariance matrix and final parameter values are used to generate 500 sample 
parameter sets. For each set, the fit fractions are calculated and recorded in histograms. Each 
histogram is fit with a single Gaussian to extract its width, which is used as a measure of the 
statistical error on the fit fraction. 



III. EXPERIMENTAL DETAILS 

The CLEO II detector is described elsewhere [^. This measurement uses the entire CLEO II 
dataset, which represents approximately 4.7 fb~^ of integrated e'^e~ luminosity at -y/i ~ 10.6 GeV. 

The D'^'s used in this analysis are required to be produced by the decay chain D^nf, 
which significantly reduces the combinatorial background. To reconstruct the D^'s, we take pairs 
of oppositely charged tracks and assign the track with the same sign as the pion from the D*~^ 
decay to be the pion from the decay. This Cabibbo-favored correlation between the signs of the 
pions eliminates the need for other particle identification techniques in this analysis. 

For tracks to be used they must be well fitted, reconstruct to within 5 cm of the interaction 
point along the beam pipe and within 5 mm perpendicular to the beam pipe (corresponding to 
about 5 standard deviations in length and more than 10 standard deviations in the width of the 
beam spot). 

We fit pairs of tracks passing these requirements to a common vertex, which is the candidate 
decay position of the meson. Each such pair of charged tracks is combined with all ir^ candi- 
dates in an event. The vr" candidates are found by combining all pairs of electromagnetic showers 
which are unmatched to charged tracks. To reduce the number of fake vr^'s from random shower 
combinations and to improve their resolution, we require that each shower have energy above 100 
MeV and be in the central region of the CLEO II detector. Furthermore, the invariant mass of the 
two photon combination is restricted to be between 128 MeV/c^and 140 MeV/c^(i.e. within about 
one standard deviation of the vr*^ mass). The two shower combination is kinematically fit to give 
the known vr*^ mass. 

Once we have a vertex with a K~ , a vr"*" and a ir^ candidate, we combine the momenta of the 
three particles to find the momentum. With the decay location and the momentum known, the 

is projected back to the beam spot. In CLEO, the beam spot has a ribbon like shape with a 
width of 700 //m, a height of 20 /im, and a length of about 2 cm. We project the candidate back 
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to the vertical position of the beam, since this dimension of the beam is most precisely known. The 
intersection of the projection and the beam position defines the production point of the D*~^. 

We refit the slow pion track to include the D*~^ production point as an additional constraint, 
providing a better measurement of its true momentum. The result of this is that the width of the 
mass difference peak, AM = M{D*~^) — M{D^), is reduced from 590 keV to 490 keV, providing a 
15% reduction in the number of background events in our final sample. We make a requirement 
that AM is between 144.9 MeV/c^and 145.9 MeV/c^. We also require that the normalized D* 
momentum, Xd* = Pd* / ij I-beam ~ ^d* ' greater than 0.6, which significantly reduces the com- 
binatorial background level and kinematically excludes the possibility that a D* candidate came 
from a decaying B meson. 

After obtaining the candidate D^'s as described above, we can plot the mass of the — > 
i^^vr+vr^ candidates as shown in Figure]^, where the fit to the mass distribution is also shown. When 
examining the Dalitz Plot, we only use the events which have 1.85GeV/c^ < M^o < 1.88GeV/c^ 
(i.e. within about one standard deviation of the known mass). 

We have chosen quite restrictive cuts on our kinematic variables (approximately one standard 
deviation on each) to minimize the effect of the background on our result. Since we are studying 
the shape of the distribution and are not trying to extract a branching ratio, the fact that this 
increases the systematic uncertainty of the overall efficiency somewhat is not an issue. 

Applying the above requirements produces 7,070 events in the Dalitz Plot. Figure ^ shows 
the distribution of this sample as a scatter plot in the chosen mass squared variables M|^_^+ and 

In order to reduce the smearing effects introduced by the detector, those combinations passing 
the above requirements are kinematically fit such that when combined, the K~ , vr"'" and vr" recon- 
struct to give the correct mass. This kinematic fit has two effects. First, the uncertainty of the 
4- momentum of the particles is reduced, giving a more precise measurement of the mass squared 
variables used to define an event's position in the Dalitz Plot. Second, the decay position in these 
variables is guaranteed to respect the kinematic boundaries of the Dalitz Plot. 



A. Background 

Turning again to Figure ^, we can see that the signal region contains a small but non-zero 
number of background events. We use the fit shown to measure the fraction of events in this 
region which are "true signal" by integrating the signal function (a double bifurcated Gaussian) 
and the background function (a line) and comparing the two. The signal fraction and its associated 
statistical error, 0.967 it 0.011, are used in the likelihood function minimized during the fitting 
procedure. 

Knowing only the amount of background is not enough if we want to correctly extract the 
amplitudes and phases of the signal component; the shape of the background in the Dalitz Plot is 
also important. There are several sidebands |jl^ that could be chosen to study the shape of this 
background using data, and a Monte Carlo study (outlined below) is used to determine which one 
is best. As will become apparent in the section on systematic errors, the overall low level of the 
background means that the final result has very little sensitivity to this choice. 

To determine which sideband will best represent the Dalitz Plot shape of the background in 
the signal region, a signal-free sample of e"'"e~ qq Monte Carlo simulated data is used (referred 
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FIG. 2. The — > K^tt^tt^ reconstructed mass distribution for all event candidates (points). 
The solid line represents a fit to the data using a double bifurcated Gaussian to represent the signal 
plus a straight line to represent the background (shaded). 



10 




M (K (GeV/c2)2 



FIG. 3. The Dalitz distribution of all 7,070 D — > K 7r"'"7r candidates in our data sam 
shown in an unbinned scatter plot. 



11 



TABLE II. Background and efficiency best fit parameters. The fitting functions are described 
in Section [V. 



Background 


Efficiency 


Bo 


1.0 (fixed) 


Eo 


(22.1 ± 1.8) X 10^5 


Bx 


-1.188 ±0.018 


Ex 


(-6.89 ±2.9) X 10-5 


By 


-0.742 ± 0.044 


Ey 


(-27.1 ±3.7) X 10-5 




U.4oo It U.uiO 


^x^ 


/iH/i-i-i (\\ N/ 1 

(^iU.4i: zt l.DJ X iU 


Bxy 


0.874 ± 0.032 


E 

^xy 


(38.2 ± 3.2) X 10-5 


By2 


0.122 ±0.047 


Ey2 


(12.4 ±2.8) X 10-5 


Bx3 


-0.052 ± 0.004 


Ex'i 


(-3.00 ± 0.27) X 10-5 


B,2y 


-0.162 ±0.010 


Ex2y 


(-7.97 ±0.73) X 10-5 


B,y2 


-0.202 ±0.014 


Exy'2 


(-12.8 ±0.94) X 10-5 


By. 


0.061 ±0.016 


EyS 


(-0.53 ±0.73) X 10-5 


B^.o 


(1.65 ± 1.70) X W-^ 




Bp+ 


(3.69 ±0.58) X 10"^ 






Bk'- 


(7.69 ± 1.95) X 10-5 







to below as the "vetoed" sample). These data are generated using a full GEANT based detector 
simulation |18|, and are processed by the same reconstruction code that is used for real data. 

This Monte Carlo sample represents the background we want to measure in data, and we use 
it as a reference in our study. The next step is to consider a number of possible sideband samples, 
and see which does the best job representing the Dalitz shape of the vetoed sample. 

Many sideband samples can be formed in the space defined by the three mass variables, AM, 
Mjjo and M^o. To choose the best one, we fit the distribution in the Dalitz Plot using an unbinned 
likelihood fit to a cubic polynomial in -M|^_^+ and M^+^o as well as non-interfering squared am- 
plitudes for the p(770), i('*(892)- and K*^ . A is formed between each Monte Carlo sideband 
sample and the reference vetoed sample to give us a measure of their relative merits. 

Based on this, the sideband which seems to best represent the vetoed sample consists of those 
events which have AM < 0.1549 GeV/c^, are in the M^o signal region, and are in the off-peak 
regions of M^o: 1.76 < M£,o(GeV/c2) < 1.80 or 1.91 < Moo(GeV/c2) < 1.95. This choice of 
sidebands along with the K^n^n^ candidates are shown in Figure ^. 

The assumption is now made that the sideband method which best represents the background 
in the Dalitz Plot when analyzing the Monte Carlo simulated data is also the best sideband method 
for use in real data. Those events from the actual data which are in the selected "best" sideband 
are fit with the cubic polynomial plus the three non-interfering resonances. The resulting best fit 
parameters are shown in Table IHj We project the fit and the background points onto the three mass 
squared variables and show the results in Figure ^, along with a two dimensional Manhattan plot 
of the fit result. We use this parameterization of the background shape in the fit to the distribution 
of events in the Dalitz Plot by including both the parameters and the covariance matrix in the final 
likelihood function (as described in Section ^ ). 
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FIG. 4. A box plot showing K vr+vr'^ candidates in the AM vs Mj^o plane. The shaded 

areas on each side of the peak are the sidebands determined to best represent the background. 
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FIG. 5. Results of the best fit to the K~ "7r~*"7r'^ sideband background sample. The fit 
function is shown as a Manhattan plot (top right), and as projections onto the three mass squared 
variables showing both fit (histogram) and data (points). 
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B. Efficiency 



Next, we determine the efficiency for detecting signal events as a function of position in the 
two dimensional Dalitz Plot. After generating 4.2 million signal Monte Carlo events with a flat 
distribution in phase space {i.e., uniform across the Dalitz Plot), these events are analyzed to find 
the number of observed events as a function of M?._ + and n. The events observed are binned 
into regions with 50 {MeY/cy on a side, and we divide the number of events observed in each bin 
by the number generated to give a measure of the efficiency for that bin. Due to the finite number 
of Monte Carlo events observed in each bin, each individual efficiency measurement has about a 
10% statistical error. Since we expect (and observe) that the efficiency is a slowly varying function 
across the Dalitz Plot, we fit the efficiency measurements with a cubic polynomial in M^_^+ and 
M^+^o and use the resulting function to parameterize the efficiency. 

As a check that the efficiency function obtained using phase space Monte Carlo is reasonable, we 
repeat the procedure described above with another 2.4 million signal Monte Carlo events generated 
with the Dalitz distribution found by E691 [||]. Again the efficiency in each bin is calculated 
and fit. Since the resulting fit agrees well with the efficiency calculated from the phase space 
distribution of points, we combine the two Monte Carlo samples and calculate the efficiency using 
the full 6.6 million events. The fit parameters for this combined fit are shown in Table Figure ^ 
shows the raw efficiency for each bin as well as the fit and the projections onto each of the three 
mass-squared axes. 



IV. FITTING PROCEDURE 

Having a parameterization for both the background and efficiency as well as knowing the fraction 
of events in the signal region which are in fact background, we can fit the data in the Dalitz Plot 
to extract the amplitudes and phases of any contributing intermediate resonances. 

To do this we use an unbinned maximum likelihood fit which minimizes the function given 

by 

■^=[ E -21n/:]+X^enalty (15) 



where 



and 



events 



^_/ F£{M],.^,,M'^,^,)\M\' ^ ^^ g(M|..^,,M^,^o) \ ^^^^ 

\ signal A/" background J 



Xpenalty = {^^^) + T,ij{Bi - Bio)Vij{Bj - Bjo) 

+ Esys Eij {Er - Eio)W,j {Ej - Ejo) (17) 

^fs^gnal = J £ [M], - ^+ , ^o) \M\^dVV (18) 

Mbackground = J l3{M^-^+ , M^+^o)dVV . (19) 
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FIG. 6. Results of the best fit to the D° K' "vr+vr detection efficiency distribution. The fit 
function is shown as a Manhattan plot (top right), and as projections onto the three mass squared 
variables of both fit (histogram) and raw efficiency (points). In each of the projections, the quantity 
plotted is the average efficiency at the given value. 
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The signal fraction F„ and its error ap (0.967 and 0.011 respectively) are determined from the 
fit to the mass spectrum shown in Figure ^, and the the parameters Bjo and Ejo describe the 
nominal background and efficiency shapes (see Table ||) via the cubic polynomial shapes 

B = Bo + S.M|_^+ + ByMl^^, + B,2(M|._^+)2 + B^yMl^^+Ml+^o + By2{M^+^of + 
BAMl-^+f + B,2y{Ml.^^)^MXo + B,y.Ml.^+{Ml+^,f + Byz{Ml+^,f + 
5^*o|A(i^"vr+7r0|;r*°)|2 + /3^|^i(^+7r0K-|p+)|2 + /?^..|^i(i^-7rV+|K*-)|2 (20) 



and 



£ = Eo + E,Ml.^+ + EyMXo + E.AMl-^+f + E^yM^^^^^M^+^o + Ey2{Ml+^,f + 

EAMl-^+f + E,2y{Ml.^+fMl+^, + E,y2Ml.^+{Ml^^,f + Eyz{Ml^^,f (21) 



In expressing this likelihood function we have made the explicit assumption that background 
events and signal events are distinct, allowing us to factor the likelihood into two components which 
do not interfere. The Xpenaity terms represent the information from the fits used to determine the 
signal fraction, the background parameterization, and the efficiency parameterization. Vij and Wij 
are the covariance matrices from the background and efficiency fits respectively. The last term is 
used only when evaluating the systematic errors due to the efficiency parameterization, hence Esys 
is set to zero during "normal" fitting. 

In addition to the likelihood, we need a measure to assess how well any given fit represents the 
data. A confidence level can be calculated directly from the likelihood function by utilizing the best 



fit parameters. This idea was described by ARGUS |20| and is a direct application of the Central 
Limit Theorem from statistics ||2^ . Assuming the candidates are truly distributed according to the 
likelihood function which gives the best fit, the average value is 



N 

i 



1 ^ r 

— ^(-21n£) « y £(-21n£)(iPP (22) 



where is the number of candidates. The variance is given by 

N 



^M = ^E(-21n/:-^)'« j C{-2\nCfdVV-^\ (23) 

i=l 

Because we have a large number of candidates distributed according to this function, the Central 
Limit Theorem tells us that the mean should follow a normal distribution. The sum of minus log 
likelihoods, which is the value minimized in the fit, has a mean of N ^ and follows a normal 
distribution with a variance of Na"^,. Thus, the minimal value will come from a normal distribution 
with mean 

<-2^\nC>=N j C{-2\nC)dVV -n (24) 

and standard deviation 

^<-2Ein/:> = V^/ ^(-21n/:)2dW - N,,^ (25) 
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where n is the number of parameters extracted from the fit. The confidence level for the fit is then 
just the area of a Gaussian with the above mean and width which hes above the value obtained 
in our fit. It is worth pointing out that this value only gives a measurement of the goodness of fit 
assuming the fit function correctly describes the true distribution. 

Having a second measure of the goodness of the fit would be extremely valuable, and an obvious 
choice is the x^- This requires the data to be binned, and furthermore that there are enough events 
in each bin that Gaussian statistics can be assumed. As we saw in Figure ^, the density of candidates 
in the —>■ K~-k^it^ Dalitz Plot varies significantly as a function of position, hence to form a 
sensible measure we will need to have bins of varying size. 

To systematically choose these bins, we start by placing a grid of small regions, 50 (MeV/c^)^ 
on a side, over the Dalitz Plot. Next, adjacent regions are combined into bins until each contains 
approximately 30 candidates. After completing this procedure, our Dalitz Plot is divided into 228 
bins of varying size, and a variable for the multinomial distribution |22,23| can be calculated as 



228 / 

X^ = -2Enan(^) (26) 



rii 



where nj is the number of events observed in bin i, and pi is the number predicted from the fit. 
For a large number of events this formulation of the becomes equivalent to the usual one |24|. 

One can naively calculate the number of degrees of freedom for the fit as the number of bins 
(r) minus the number of fit parameters (k) minus one, as would be correct for a binned maximum 
likelihood fit. However, since we are minimizing the unbinned likelihood function, our "x^" variable 
does not asymptotically follow a x^ distribution |24], but it is bounded by a x^ variable with (r — 1) 



degrees of freedom and a x^ variable with {r — k — 1) degrees of freedom. Because it is bounded 
by two x^ variables, it should be a useful statistic for comparing the relative goodness of fits. In 
what follows, we use both the x^ ^-nd the confidence level described above as our "goodness of fit" 
measures to determine which of the many possible sets of intermediate resonances are preferred. 

Before analyzing the — > K~t^^-k^ data, we performed many checks of both the fitting and fit 
evaluation procedures. One of these was a double-blind study in which several Monte Carlo samples 
containing — > K^-k^tt^ decays generated with "secret" mixtures of intermediate resonances were 
analyzed. In each case, our fitting and evaluation procedure identified the correct set of resonances, 
and recovered their amplitudes and phases within statistical errors. The resulting amplitudes and 



phases for one of the fits is shown in Table [II . 



V. FITTING THE DATA 

Armed with the tools described in the previous section, we are ready to fit the data distribution 
shown in Figure 0. Previous experiments have observed three intermediate resonances in 
K~'K^-n'^ decays: /O"*", K* and K*~ , hence we begin by considering only these in addition to a 



non-resonant component. The resulting fit parameters are given in Table IV. 

Figure ^ shows the projections of both the fit and the data onto the three mass squared variables, 
as well as a two dimensional Manhattan plot of the final fit function. Even a quick glance suggests 
that the data is not well represented by this function, and the large value of x^ as well as the 
zero confidence level confirm this observation. These parameters are useful for comparison with 
previous experiments, however, which reported observation of these three resonances with much less 
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TABLE III. A comparison between input Monte Carlo parameters and the results from a 
subsequent fit to the Dalitz Plot using the techniques described in Section IV. Note that the input 
amplitudes and phases are completely fictitious. 



Resonance 




Generated 




Measured | 




Amplitude 


Phase (degrees) 


Amplitude 


Phase (degrees) 




1.0 


45 


1.03 ± 0.02 


47 ± 1 




1.0 





1.0 (fixed) 


(fixed) 




1.0 


-115 


1.03 ±0.02 


-113 ±2 


i^o*(1430)- 


0.5 


-115 


0.54 ±0.05 


-107 ±6 


Non resonant 


1.0 


-90 


1.08 ±0.05 


273 ±3 



TABLE IV. Results of the best fit to the data with only p^, K , K* , and non-resonant 
components included. 



0>nr 


1.70 ±0.07 


ap+ 


1.00 (fixed) 


ttK*- 


0.378 ± 0.008 




0.422 ± 0.009 


4>nr 


59.7° ± 2.0° 


<l>p+ 


0° (fixed) 


4'K*- 


166.7 ±2.0° 




-7.8° ±2.2° 


-21n£ 


7070 


Conf. Level. 


0.0% 




650 
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TABLE V. A comparison of the fit fractions obtained with our "three resonance" fit and those 
reported by previous experiments. The errors shown are statistical only. Note that although 
the data is not well fit by this model, the results are consistent with those reported by previous 
experiments. 



Decay Mode CLEO II (3 Resonance) E687 Mark III E691 

K-p+ 0.834 ± 0.007 0.765 ± 0.041 0.81 ± 0.03 0.647 ± 0.039 

K*-TT+ 0.129 ±0.006 0.148 ±0.028 0.12 ±0.02 0.084 ±0.011 

K*^TT^ 0.157 ±0.007 0.165 ±0.031 0.13 ±0.02 0.142 ±0.018 

Non resonant 0.074 ± 0.006 0.101 ±0.033 0.09 ± 0.02 0.036 



TABLE VI. A comparison of the phases (in degrees) obtained with our "three resonance" fit and 
those reported by previous experiments. The errors shown are statistical only. In the "Rotated" 
column we have shifted the reported phases such that the p has a phase of 0° in order to ease 
comparison with the other results. Note that although the data is not well fit by this model, the 
results are consistent with those reported by previous experiments. 



Decay Mode CLEO II (3 Resonance) E687 Mark III E691 (Rotated) 

K-p+ (fixed) (fixed) (fixed) 0± 7 

K*-TT+ 166.7 ±2.0 162 ±10 154 ±11 -152 ±9 

-7.8 ±2.2 -2 ±12 7±7 127 ± 9 

Non-resonant 59.7 ± 2.0 -122 ± 10 52 ± 9 -40 (fixed) 



statistics. We show the comparison in Tables ^ and VI and see good agreement. Unfortunately, we 
can only compare the results for the phases to E687 since the other experiments do not give their 
choice of particle ordering or potential complex constants in their choice for A{ABC\r). Although 
the phases match for the three resonant components, the non-resonant phase seems to be off by 
180°. This observation is consistent with comments that E687 had an unreported negative sign in 
their vector amplitude |25]. 

Since we have at least a factor of ten more statistics for this analysis, one should not be surprised 
that more resonances are needed to accurately represent the data. The question now becomes how 
best to determine which additional resonances to include. We have tried two procedures: a) adding 
all possible resonances and subsequently removing those which do not contribute significantly, and 
b) adding new resonances one at a time and choosing the best additional one at each iteration, 
stopping when no additional resonances contribute significantly. Both of these methods lead us to 
the same results, hence only the first one is described below. 

We begin by fitting the Dalitz Plot with all known resonances which can possibly contribute 
to this decay, as listed in Table VII [2^ . The results of this fit are shown in the "All Resonances" 
column of Table VIII, and in Figure ^. There are five resonances which have fit fractions that are 
less than one standard deviation away from zero: ^3(1780)°, ir|(1780)~, Z*(1410)°, /r*(1410)- 
and i^*(1680)*'. Two other resonances, -fC2(1430)*^ and ii'|(1430)~, have fit fractions close to zero. 
When the first five resonances are removed and the fit repeated, the fit fractions of these last two 
resonances do become consistent with zero, and hence are also removed. 



Notice that in the "All Resonances" column of Table VIII there are two heavy p mesons 
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3111100-016 




FIG. 7. The results of fitting the K' 7r"'"7r'^ data with only K , K* , and non 

resonant components included. The efficiency corrected fit is shown as a Manhattan plot (top 
right), and as projections onto the three mass squared variables of both fit (histogram) and data 
(points). The dashed line shows the level of the background. 
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FIG. 8. The results of fitting the K' 'tt~^it data with all known resonances likely to 
contribute (see Table VII). The efficiency corrected fit is shown as a Manhattan plot (top right), 
and as projections onto the three mass squared variables of both fit (histogram) and data (points). 
The dashed line shows the level of the background. 
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TABLE VII. The resonances considered when fittmg the K vt+tt'^ Dahtz Plot, along 

with the masses and widths used when evaluating the matrix element. 





Parameters 


Resonance 


J 


Mass (LreV /c j 


Tier- Jj-U /r^^\T /^2\ 

Wiatn (LreV /c j 


/)(770)+ 


1 


n '7'7n 1 r\ r\ni 

0.770 ± 0.001 


0.1507 ± 0.0011 


(892)0 


1 


O-oyol ± O.OOOo 


0.0505 ± 0.0006 


i^*(892)- 


1 


0.8915 ± O.OOOo 


0.050 ± 0.001 


i^*(1410)" 


1 


1.414 ± 0.015 


0.2c!2 ± 0.021 


1410 


1 


1.414 ± (J.Ulo 


(J.zJz ± U.Uzl 


i^o*(1430)~ 


0+ 


J-.^iJ-Zi \J.VJVJU 




(1430)0 


0+ 


1.412 ±0.006 


0.294 ± 0.023 


i^l (1430)- 


2+ 


1.425 ± 0.002 


0.098 ± 0.003 


(1430)0 


2+ 


1.432 ±0.001 


0.109 ± 0.005 


/)(1450)+ 


1" 


1.465 ± 0.025 


0.310 ±0.060 


/9(1700)+ 


1- 


1.700 ± 0.020 


0.240 ± 0.060 


e:*(1680)- 

(1680)0 


1" 


1.717 ±0.027 


0.322 ±0.110 


l- 


1.717 ±0.027 


0.322 ±0.110 


(1780)0 


3" 


1.776 ± 0.007 


0.159 ±0.021 


i^3*(1780)- 


3" 


1.776 ± 0.007 


0.159 ±0.021 



(/9(1450)^ and p{1700)~^) which have surprisingly large fit fractions. Both have masses which 
place their peak outside the Dalitz Plot, but both are wide enough (310 ± 60 MeV/c^ and 240 ± 60 
MeV/c^ respectively |^6|) that their tails extend well into the region of interest, making it difficult 
to distinguish between them. Since the fitted phases of these /o's are very close to being 180° apart, 
their large fit fractions are assumed to be an artifact of the fit's inability to tell them apart. Sup- 
porting this is the additional fact that when both resonances are combined, their net contribution 
to the fit fraction is much smaller, (9 ± 2)%. Since the inclusion of both p resonances is probably 
a misrepresentation of the contents of the Dalitz Plot, only one of these is included in all following 
fits. We choose the one which gives the best ™d goodness of fit, the ^(1700)"*", and consider the 
/9(1450)^ only when evaluating our systematic errors. 

After the seven resonances consistent with zero fit fraction are removed along with the /9(1450)~'~ 
(as discussed above), seven resonances remain in addition to the non-resonant component: /?(770)^, 
i^*(892)-, (892)0, p(i7oo)+, ;^o(1430)0, Ko(1430)^, and i^*(1680)-. Figure | shows the result 
of fitting the Dalitz Plot with these components. The fit fractions and phases are shown in the 
"Final Resonances" column of Table VIII , and the full set of parameters extracted from this fit is 
shown in Table [X. 

As a curious side note, if a single vector {K~tt^) resonance with a floating mass and width is 
added in place of the four new "standard resonances" discussed above, a good fit can obtained. 
Unfortunately, while this new resonance has a reasonable mass of 1.406 GeV/c^, it prefers a negative 
width of r = —0.25 GeV/c^ which does not seem to represent the underlying dynamics we are trying 
to measure. It is possible that the desire for this resonance is an indication of an inaccuracy of 
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TABLE VIII. The parameters from tlie fits to tlie D K^tt^tt Dalitz Plot with ah reso- 
nances included ("All Resonances" column), and after we remove resonances consistent with zero 
fit fraction ("Final Resonances" column). The /9(1450)+ and p(1700)+ contributions are discussed 
in the text. 





All Resonances 


Final Resonances 




J. IICXOC \V_iC&jiCCO, 




Fit Frartinn (%) 


FVirmp (Hpp"rpp<*l 




Fit Fraction (%\ 


^3(1780)° 


263 ± 16 




0.3 ±7.5 










86± 12 




0.5 ±2.9 








i^* (1680)0 


175 ± 25 




0.4 ±0.5 








X*(1680)- 


67± 19 




1.0 ±0.5 


103 ±8 




1.3 ±0.3 


p(1700)+ 


149 ±8 




75± 18 


171 ±6 




5.7 ±0.8 


,9(1450)+ 


-45 ± 10 




34± 11 








Non Res. 


30 ±5 




9.1 ± 1.3 


31 ±4 




7.5 ±0.9 


ir*(1410)0 


279 ± 52 




0.1 ±0.2 








^2 (1430)0 


148 ± 13 




0.3 ±0.14 








(1430)0 


168 ±5 




8.0 ± 1.3 


166 ±5 




4.1 ± 0.6 


K*(1410)- 


152 ± 31 




0.2 ±0.2 








i^2(1430)- 


339 ± 21 




0.12 ±0.08 








i^o(1430)- 


42 ±6 




5.6 ± 1.1 


55.5 ± 5.8 




3.3 ±0.6 


K*{892)- 


159 ±2.6 




12.8 ± 1.8 


163 ±2.3 




16.1 ±0.7 


p(770)+ 


(fixed) 




74 ±4 


(fixed) 




78.7 ±2.0 


i^* (892)0 


2.8 ±3.2 




11.3 ± 1.5 


-0.2 ±3.3 




12.7 ±0.9 






203 






257 




-2 In/: 




6490 






6570 




C.L. 




91.3 






94.9 
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TABLE IX. Summary of our best fit to the data with the final set of eight components included. 



Sigiiiil Piiramel 




Background Paramel crs 




1.75 ±0.12 






ap+ 


1.00 (fixed) 


Bo 


1.0 ±0.0 




0.44 ±0.01 


Bx 


-1.206 ±0.001 




0.39 ±0.01 


By 


-0.74 ± 0.23 




0.77 ± 0.08 




0.468 ± 0.001 


"k'o (1430)0 


0.85 ± 0.06 


Bxy 


0.842 ± 0.008 


ap(1700)+ 


2.50 ±0.19 


By2 


0.168 ±0.001 




2.50 ±0.3 


B^3 


-0.055 ±0.001 


4>NR 


31.2" ±4.3° 


B 2 

±J^^y 


-0.16 ±0.06 


<f>p+ 


0" (fixed) 


Bx,. 


-0.188 ± 0.001 


4>K*- 


163 ± 2.3° 


By. 


0.077 ± 0.001 




-0.2° ±3.3° 


B-*o 


(3.4 ±0.1) X 10-5 


4'Kq(IAZQ)- 


55.5° ± 5.8° 


Bp 


(4.27 ±0.05) X 10"^ 


'%o (1430)0 


166° ± 5° 


Bk*- 


(9.64 ±0.01) X 10-5 


0p(17OO)+ 


171° ± 6° 






^^♦(leso)- 


103° ± 8° 






Signal Fraction 


0.968 ± 0.007 






-21n£ 


6570 






Conf. Level. 


94.9% 








257 
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FIG. 9. The results of fitting the K^tt^tt^ data with the "final set" of components. 

The efficiency corrected fit is shown as a Manhattan plot (top right), and as projections onto the 
three mass squared variables of both fit (histogram) and data (points) . The dashed line shows the 
level of the background. 
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the formalism used for the resonance shapes, or an indication that multiple resonances are needed 
(as we have assumed). We note that the optimum set of seven resonances used above, all of which 
have positive widths, provide a fit which has a lower than the inclusion of this single unphysical 
state. 

Other experiments have reported evidence of a light scalar (7r"'"7r~) resonance, the a, in 
ir^ir^ir^ decays |2^, as well as evidence of a scalar {K^ir'^) resonance, the k, in K'-k^tt^ 
decays ||2^. Since a significant fit fraction for L)+ kvt"'" has been reported by these authors, we 
have searched for a scalar k — > K~7:~^ resonance in the — > ktt^ channel, fixing the mass and 
width of the k to the values reported in |^, (0.815 GeV/c^ and 0.560 GeV/c^ respectively). We 
find a fit fraction consistent with zero (0.4 it 0.3%), and see no improvement in the confidence level 
of the fit with this additional resonance included. We have also allowed the mass and width of the 
K to float in the fit, and again see no significant contribution. 

Lastly, since this analysis considers only mesons produced from a decaying D*~^ in the mode 
D*^ D^irf, we have the ability to divide our data into separate and samples by simply 
considering the sign of the vr^ from the D*^ decay. The Dalitz Plots of these samples can then 
be fitted separately and compared in a search for CP violation. We have fitted these samples 
with the same set of resonances described above, and the results are shown in Table When 
performing these fits, the efficiency functions were found separately for the and samples, 
however a common background shape was assumed. Forming a simple between the two sets of 
fit parameters we find Xcp = 16.2 for 14 degrees of freedom. 

We calculate an integrated CP asymmetry across the Dalitz Plot by evaluating 

and obtain Acp = —0.031 it 0.086, consistent with zero. Note that this number is not dependent 
on the number of and candidates in our data sample, but rather on the shapes of these 
distributions in the respective Dalitz Plots. 



VI. SYSTEMATIC UNCERTAINTIES 

After finding the best fit to the data, we must attempt to estimate the systematic uncertainties 
on the fit parameters. There are several possible sources: the background, the efficiency, biases 
due to experimental resolution, and the modeling of the decay. These contributions are discussed 



in order, and the final systematic errors are shown in Table Xl, where experimental and model 
dependent sources of systematic uncertainty are summarized in detail. 

The background was modeled by the choice of sideband sample that gave the best parameter- 
ization of the vetoed data sample from Monte Carlo. Furthermore, the background parameters 
were allowed to float in our fits to the data, constrained only by the covariance matrix from the 
fit that determined the nominal background function. To search for any systematic effects due the 
background parameterization, the fitting procedure was repeated for a number of different sideband 
choices. Because the background fraction our sample is a mere 3.3%, or about 230 out of the 7070 
events in the Dalitz plot, these changes have a minimal effect on the fit parameters. We use the 
RMS spread of these results as our estimate of the systematic error due to our choice of background 



parameterization. These values are shown in the "Bkgnd" column of Table XI 
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TABLE X. Fit results when the and samples are considered separately. 





1)0 


Sample 


Sample 


Component 


Amplitude 




Fnase (degrees) 


Amplitude 




Fnase (degrees) 


p(770)+ 


1.0 ±0.0 




(nxed) 


1.0 ± 0.0 




(nxed) 


K*{892)- 


0.433 ± 0.034 




168.9 ± 3.3 


0.442 ± 0.015 




157.8 ± 3.4 


(892)0 


0.391 ± 0.026 




1.3 ±3.7 


0.410 ± 0.022 




-4.9 ±4.9 


/9(1700)+ 


2.590 ±0.538 




175.0 ±7.5 


2.720 ±0.272 




163.9 ±7.6 


Zo(1430)° 


0.989 ±0.124 




173.9 ±8.2 


0.774 ± 0.089 




159.3 ±8.1 


i^o(1430)- 


0.701 ±0.211 




59.0 ± 10.0 


0.917 ±0.117 




55.0 ±7.1 


i^*(1680)- 


2.567 ± 1.540 




107.4 ± 69.2 


2.060 ± 0.423 




106.4 ± 13.5 


Non Res. 


1.840 ±0.146 




39.9 ±7.9 


1.780 ±0.160 




21.3 ±6.0 






227 






233 




-21n£ 




3237 






3302 




C.L.(%) 




93.1 






80.7 





To obtain the efficiency across the Dalitz Plot, signal Monte Carlo events were fit to a cubic 
polynomial. As a check, we have allowed this polynomial to float in our fit (as was done with 
the background) subject to a constraint from its covariance matrix in the likelihood function 
(i.e. setting Esys = 1 in Equation 17). If the efficiency is not well modeled by a cubic polynomial, 
there could still be an effect that this check would fail to find. To search for this we tried a local 
smoothing algorithm rather than the global polynomial fit. The efficiency was smoothed by fitting 
either nine or twenty-five neighbors around each bin with a local plane. Each bin's efficiency value 
was then replaced by the height of this plane interpolated to its center. As a final check, we 
used the raw measurements of the efficiency in each bin of our fits. We conclude that the effects of 
parameterization of the efficiency function over the Dalitz Plot is not a significant source of concern 
as most of the fit parameters vary by less than their one sigma error bars in the above checks. 

Since we make no requirement on the momentum of the charged tracks, one might worry that 
low momentum tracks may be poorly measured and could affect the Dalitz Plot distribution in a 
way not well modeled by our Monte Carlo. To search for such a momentum dependent effect, we 
fit the data with the additional requirement that all tracks have a momentum above 350 MeV/c. 

The cuts used to obtain our signal determine the structure of our efficiency. To assess how well 
the Monte Carlo reproduces the data distributions, we varied the cuts used in the analysis and fit 
the resulting Dalitz distributions. Each cut was relaxed in turn. The cuts on the masses, Mjjo, 
AM and M^o, were opened to double the size of the signal region. The minimum energy on the 
photons was relaxed to 90 MeV, and the requirement on X^* was loosened to 0.5. 

The RMS variation in the fit parameters from each of the tests described above was taken as 
our estimate of the systematic uncertainty on the efficiency. These values are shown in the "Eff" 
column of Table KI. 

A final contribution to the experimental systematic error, presented in column "Resol" of Ta- 
ble Xl, is due to the finite resolution of the Dalitz Plot variables. As a check, we have included 
the effects of smearing when fitting the data. This was done by measuring the resolution as a 
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TABLE XI. A summary of the systematic errors on each fit parameter. The first two columns 
show the results from the best fit and the associated statistical errors. The next four (three) 
columns summarize the systematic uncertainties due to experimental (modeling) sources respec- 
tively. Details are provided in the text. 





From Fit 


Experimental Errors 


Modeling Errors 


Parameter 


Value 


Stat, J^iVv 


Bkgnd 


Eff 


R psol 


Total 


Sh rlDP 


Add 


Total 


(892)0 Frac (%) 




Q1 


n 1 7 






U.(JU 


1 49 


n "^s 

u.oo 


1 47 


Z* (892)0 pYias^ (deg.) 




o. Z(0 


1 nfi 

l.UU 


1 fi9 


1 fizL 


9 9n 


R QQ 


n fi7 


7 09 


/)(770)+ Fit Frac (%) 


78 7fi 




n ^9 


1 1 n 


u.oo 


1 .oo 


4 4n 


x.oo 


4 fin 


i^*(892)- Fit Frac (%) 


16.11 


0.69 


0.47 


0.53 


0.18 


0.73 


+2.58 
-0.48 


0.59 


+2.65 
-0.76 


i^*(892)- Phase (deg.) 


163.40 


2.32 


0.94 


2.62 


1.30 


3.08 


4.20 


1.09 


4.34 


Ko{U30)~ Fit Frac (%) 


3.32 


0.64 


0.13 


0.60 


0.40 


0.73 


1.16 


0.40 


1.23 


Ko(1430)- Phase (deg.) 


55.52 


5.76 


1.20 


2.76 


1.31 


3.28 


-12.8 
+3.1 


2.85 


+4.2 

-13.1 


(1430)0 Fit Frac (%) 


4.05 


0.61 


0.15 


0.66 


0.24 


0.72 


+3.04 
-0.24 


0.39 


+3.06 
-0.46 


Ko(1430)0 Phase (deg.) 


165.90 


5.23 


2.39 


3.83 


0.70 


4.57 


11.4 


3.20 


11.8 


p(1700)+ Fit Frac (%) 


5.65 


0.76 


0.20 


0.43 


0.50 


0.68 


5.71 


0.59 


5.74 


/9(1700)+ Phase (deg.) 


170.50 


6.07 


1.90 


3.90 


1.50 


4.59 


-54.7 
+3.3 


5.17 


+6.1 
-54.9 


i^*(1680)- Fit Frac (%) 


1.33 


0.33 


0.07 


0.32 


0.11 


0.34 


0.17 


0.32 


0.36 


K*(1680)- Phase (deg.) 


103.20 


7.90 


3.71 


5.91 


2.00 


7.26 


9.21 


9.89 


13.5 


Non Res Fit Frac (%) 


7.50 


0.95 


0.35 


0.42 


0.05 


0.55 


+5.54 
-0.79 


0.41 


+5.56 
-0.89 


Non Res Phase (deg.) 


31.20 


4.28 


1.28 


5.08 


1.70 


5.51 


-14.4 
+3.5 


1.19 


+3.7 
-14.4 



function of position across the Dalitz plot and numerically convoluting this with the amplitude at 
each point when performing the fit. Again, the parameters vary by less than the statistical errors 
on the nominal best fit, and their variation from the nominal values is taken as an estimate of the 
systematic uncertainty. 

The above three systematic error categories (background, efficiency and resolution) are sum- 
marized in Table KI. They are combined in quadrature to give the total experimental uncertainty, 
which is shown in the "Total" column under "Experiment" . 

Modeling systematic errors can arise from our choice of resonances and the uncertainty in their 
shapes. In Section |ll| we motivated our choice of parameterization of the intermediate resonances; 
however, other groups have used different functional forms in their fits |§J^. We varied these shapes 
to study any systematic effects resulting from our choice. We examine three variations: (i) the 
Zemach formalism which enforces the transversality of the mesons by using rather than 

in the denominator of the spin sums, (ii) a simple cosine distribution for the spin sum and (iii) a 
non-relativistic rather than relativistic Breit-Wigner in the propagator. Further consideration was 
also given to the radial parameters used in the form factors, which were varied between GeV^^ 
and 3 GeV~^ for the intermediate resonances and between GeV"^ and 10 GeV"^ for the 
meson. The masses and widths of the intermediate resonances were allowed to vary within the 
known errors [26|. The non-resonant contribution was described in our fits by a constant term, but 
as a check we also modeled it by a linear function or a shape given by the spin structure without 
the Breit-Wigner amplitudes [31|. 
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The above tests were used to explore the systematic dependance of the fit parameters on the 
way the physics was modeled. The variations using a simple cosine distribution in place of the 
spin sum and using a spin structured rather than constant non-resonant component resulted in fits 
with significantly worsened (368 and 322 respectively), and are not considered when assigning a 
systematic error as the data suggests these forms could not be correct. We take the largest of the 
remaining variations as the systematic error due to our choice of modeling shapes, and the results 
are shown in the "Shape" column of Table KI. 

The final systematic check is on our choice of which resonances to include. For example, there 
is only a slight preference for the p(1700)''' over the p(1450)''' based on the goodness of fit. To 
account for this uncertainty, both fits were performed and the variation of the parameters were 
noted. Fits were also performed which included additional resonances from Table VII . The RMS 
variation in the fit parameters from the above checks is presented in the "Add" column of Table Kl . 

We also considered the effects of removing resonances, and two of these studies deserve further 
comment. The first is the removal of the K*{1680)~. We considered this because the final fit 
fraction for this resonance is a rather small 1.3 it 0.3%. When the i^*(1680)~ is removed the 
increases from 257 to 316 indicating that this resonance should remain. The parameters for this 
fit are shown in the "Removed /C*(1680)~" column of Table XII. For comparison, when the other 
"new" resonances, Kq{1430)^ , i^o(1430)~, and p(1700)"'", are removed, the increases to 379, 
348, and 381 respectively. The second case which deserves special attention is the removal of the 
non-resonant 

component. Some theoretical models, such as chiral perturbation theory ||3^, prefer a small 
non-resonant component, suggesting it proceeds only by the coherent sum of two body decays. 
When this test is performed on our data, the resulting jumps to 411, suggesting that a non 
resonant component is indeed present. The parameters for this fit are shown in the "Removed Non 
Resonant" column of Table XD . 

Since removal of any of the fit components causes a significant increase in the of fit, these 
variations were not included in the modeling systematic error. To obtain the final model dependent 
systematic error we add the "Shape" and "Add" columns of Table XI in quadrature to obtain the 
result shown in the "Total" column under "Model" . 



VII. SUMMARY OF RESULTS 

We have fit the distribution of data in the K^ir^ir^ Dalitz Plot obtained with the CLEO II 

experiment to a coherent sum of seven intermediate resonances plus a non-resonant component. 
All resonances are either scalar or vector; no significant tensor contribution was found. The non- 
resonant contribution is significant, and cannot be removed without seriously compromising the 
quality of the fit. We see no evidence of a scalar k K~7r~^ resonance in the mass range recently 
reported by other groups. 



The final fit fraction and phase for each component is given in Table XIII, These fit fractions, 
multiplied by the world average ^ K- '7r~^7:^ branching ratio of (13.9 it 0.9)% [p7|, yield the 



partial branching fractions shown in Table XIV. The error on the world average branching ratio 



is incorporated by adding it in quadrature with the experimental systematic errors on the fit 
fractions to give the experimental systematic error on the partial branching fractions. Note that 
due to interference the fit fractions do not add to unity, and consequently the partial branching 
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TABLE XII. Fit results after removal of the either the i^*(1680) resonance or the 
non-resonant component. See Section VI for discussion. 





Removed K*{1680)- 


Removed Non Resonant 


Component 


Phase (degrees^ 




r It rraciion {/o) 


Phase (degrees^ 




riz ri action {/o) 


p(770)+ 


U (hxedj 




8U.8 ± 8.5 


U(hxeaj 




I i.o ± i.8 


i^*(892)- 


157 ±6.7 




13.8 ± 1.0 


161 ±2.2 




18.2 ±0.7 


:^(892)0 


-4.7 ±5.7 




14.5 ± 1.3 


-1.5 ±2.8 




10.7 ±0.8 


/9(1700)+ 


161 ± 20 




6.7 ±0.8 


161 ± 5 




5.4 ±0.8 


Zo(1430)° 


164 ±9 




4.4 ±0.5 


194 ±9 




1.0 ±0.3 


i^o(1430)- 


47.8 ±3.6 




4.5 ±0.7 


11 ±4 




5.2 ±0.7 


K*{W80)- 







0.0 


90 ±5 




1.9 ±0.5 


Non Res. 


37 ±6 




7.7 ± 2.6 







0.0 






316 






411 




-21n£ 




6653 






6798 




C.L.(%) 




98.5 






0.7 





fractions do not sum to the world average. 

By separately fitting the — > K~7T~^7r^ and K^tt^tt^ Dalitz Plots, we have calculated 

the integrated CP asymmetry across the Dalitz Plot to be Acp = —0.031 ± 0.086. 
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ration, A.H. Mahmood thanks the Texas Advanced Research Program, F. Blanc thanks the Swiss 
National Science Foundation, and E. von Toerne thanks the Alexander von Humboldt Stiftung for 
support. This work was supported by the National Science Foundation, the U.S. Department of 
Energy, and the Natural Sciences and Engineering Research Council of Canada. 



TABLE XIII. Final fit results. The errors shown are statistical, experimental systematic, and 



modeling systematic respectively, as discussed in Section VI and summarized in Table XI 



Mode 


Fit Fraction 


Phase (degrees) 


p{770)+K- 

Z*(892)07r0 

p(1700)+i^~ 

Zo(1430)07rO 

K^(1430)-^+ 

ivr*(1680)-7r+ 

Non Resonant 


0.788 ± 0.019 ± 0.013 ± 0.046 

0.161 ±0.007 ±0.007l[!;[3o8 
0.127 ± 0.009 ± 0.005 ± 0.015 
0.057 ± 0.008 ± 0.007 ± 0.006 

0.041 ± 0.006 ± 0.007l[!:J]|j^ 
0.033 ± 0.006 ± 0.007 ± 0.012 
0.013 ± 0.003 ± 0.003 ± 0.003 

0.075 ± 0.009 ± 0.006t[!:[]^^ 


0.0 (fixed) 
163 ±2.3 ±3.1 ±4.3 
-0.2 ±3.3 ±2.2 ±7.0 
171 ± 6 ± bllti 
166 ± 5 ± 4.6 ± 12 
55.5 ± 5.8 ± 3.31^3^ 
103±8±7± 14 
31 ± 4 ± 5.5tl\ 



31 



TABLE XIV. Partial branching fractions calculated by combining onr fit fractions with the 
previously measured K~7r~^Tr^ branching ratio as described in the text. The errors shown 

are statistical, experimental systematic, and modeling systematic respectively. 



Mode 



Partial Branching Fraction 



B{D^ 

B{D^ 
B{D^ 
B{D^ 
B(D^ 



p{770)+K-) X B{p{770)+ - 
K*(892)-7r+) X B{K*{892) 
X*(892)07r0) 
p(1700)+i^- 
Ko(1430)07rO 



vr ' 



; B(K* {892f ^ K-Tr+) 
X B(p(1700)+ ^ 7r+7r°) 

X B(Kl{imf 



IT ' 



Ko*(1430)-7^^ 
K*(1680)-7r^ 



X 5 (K* (1430)- ^ K' 
X B(K*{168Q)- K- 



IT 



K 7r"^7r ) Non Resonant 



0.109 ± 0.003 ± 0.007 ± 0.006 



0.022 ± 0.001 ± 0.002: 



-0.004 
-0.001 



0.018 ± 0.001 ± 0.001 ± 0.002 
0.008 ± 0.001 ± 0.001 ± 0.001 



0.006 ± 0.001 ± 0.001 



-1-0.004 
-0.001 



0.005 ± 0.001 ± 0.001 ± 0.002 
0.0018 ± 0.0004 ± 0.0004 ± 0.0004 

0.010 ± 0.001 ± o.ooil[j:[jg^ 
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